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ABSTRACT 

We show that the optical counterpart of the X-ray source CX 1 in M4 is 
a ~20th magnitude star, located in the color-magnitude diagram on (or very 
close to) the main sequence of the cluster, and exhibiting sinusoidal variations 
of the flux. We find the X-ray flux to be also periodically variable, with X- 
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j^ ■ ray and optical minima coinciding. Stability of the optical light curve, lack of 

UV-excess, and unrealistic mean density resulting from period-density relation 
for semidetached systems, speak against the original identification of CX 1 as a 
cataclysmic variable. We argue that the X-ray active component of this system 
^ . is a neutron star (probably a millisecond pulsar). 

vn 

(T^ \ Subject headings: binaries: close - globular clusters: general - globular clusters: 

^ \ individual (M4) - cataclysmic variables - X-rays: stars 



^ i 1. Introduction 



Globular clusters (GCs) contain a large population of exotic binary stars with degen- 
erate components, and many of these are X-ray sources with various degrees of activity. 
Numerous such sources have been detected by the Chandra space observatory in over 80 
distinct clusters, pro viding a valu able tracer population for large-scale numerical simulations 



of cluster evolution (jPooleyl 120101 ) . 



*Based on observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hub- 
ble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), 
the Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data 
Centre (CADC/NRC/CSA). 
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One of the clusters surveyed by Chandra is M4 (NGC 6121). Eight of the X-ray sources 
found in its core region were classified as chromospherically active binaries, one as a cat- 
aclysmic variable (CV), and one was identified with the millisecond pulsar PSR1620-26 
(JBassa et al.l 12004 hereafter BPH). The CV-candidate, referred to as CX 1 in BPH, is the 
brightest among them (Ivo.5-2.5keV = 8.3 x 10^^ erg s~^). BPH identified it with a bright 
V555 = 17.37 star, located on the main sequence of M4 just below the turnoff. They argued 
that its X-ray luminosity is too high for a magnetically active binary composed of two main- 
sequence stars (i.e., a BY Dra system), while its optical luminosity is too low for an active 
binary with a subgiant (i.e., an RS CVn system). Their identification of CX 1 as a CV was 
additionally supported by its rather hard X-ray spectrum (i^o.5-i.5kev/-^i.5-6.okeV = 1-22), 
and a high ratio of X-ray and optical fluxes. 

BPH noted that the optical centroid of the bright star is offset by 2a from the Chandra 
position of CX 1, so the identification might be questionable. If that were the case, the 
actual counterpart of CX 1 could be much fainter, possibly lost in t he glare of the b r ighter 
star. Indeed, the examination of the M4-photometry published by I Anderson et al.l (120081 ) 
showed that the object identified by BPH was a blend of two stars. In the present Letter 
we demonstrate that the X-ray activity is associated with the much fainter and redder 
component of the blend. We argue that CX 1 is a millisecond pulsar or qLMXB rather than 
a cataclysmic variable. 



2. Observations 



2.1. Ground-based 



M4 was observed in nine seasons betw een 1995 and 2009 as one of the objects included 
in the CASE project (JKaluzny et al.ll2005l ). Altogether, 1683 y-band images were collected 
with the 2.5-m du Pont telescope at Las Campanas Observatory, using the same CCD camera 
and the same filter. Exposure times ranged from 20s to 120s, with a median value of 45s. The 
median seeing was 0.98". The photometric analysis was based on a combination of image sub- 
traction and profile ph otometry techniques. We used scripts emp loying the Daophot/A Ustar 
package (Stetsonlll987l) and a modified version of the ISIS code (JAlard &: Luptonlll998l ): see 
Kaluzny et al.l ( 20101 ) for details. Using ISIS, one obtains a light curve in the form of differ- 
ential counts. This allows for the photometry of variable objects even in those cases when 
they are tightly blended. Conversion of the light curve from differential counts to stellar 
magnitudes requires determination of the pedestal flux for a target star on the reference 
image (see Sect. 12.21) . Further details concerning the observations and the calibration of the 
photometry will be given in a forthcoming paper (Kaluzny et al. 2012; in preparation). 
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The optical counterpart of CX 1 proposed by BPH turned out to be a low-amplitude 
periodic variable with a stable light curve, suggestive of an ellipso idal and /or irradiation 
effect . The analysis of variance, implemented in the Tatry code (ISchwarzenberg-Czerny 
19961 ). yielded a period Pi = 0.2628216(1) d. The linear ephemeris 



HJD,^i^ = 2, 454, 656.7394(16) + E x P^ 



(1) 



provides a good fit to all data - we found no evidence for any change of period, amplitude 
or average magnitude from one observational season to another. There is, however, an 
ambiguity concerning the period: a value P2 = 0.5256432(1) d, i.e. two times larger than 
Pi, is also perfectly acceptable (see Fig. [T]). We address this issue in Sects. 12.21 and 151 
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Fig. 1. — Left column: optical light curves of the star #5274 phased with the ephemeris 
([T]). a: ground-based, \^-band (mean error of a single observation is 0.20 mag); b: the same 
data binned (50 bins per period); c: archival HST data, F435W-band, ACS/ VEGA system; 
d: color index (F435W - V^). Right column: the same data phased with P2 used in eq. ([T]) 
instead of Pi . 
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Fig. 2. — Location of stars #5273 (triangle) and 7^5274 (point with errorbars) on the CMD 
of M4. All magnitudes are given in the ACS/VEGA system. The F658N filter is centered 
on the Ha line; the F606W magnitudes are interpolated to the phases at which frames in 
filters F814W, F435W and F658N were tak en (0.666, 0.513 and 0.073, respectively). The 
F606W-F814W CMD is taken directly from lAnderson et al.l (120081 ) . and the remaining two 
CMDs are extracted from subfields of the images mentioned in Sect. 12.21 



2.2. HST 



The optical counterpart of CX 1 is visible not only in HST/WFPC2 images analysed 
by BPH, but also in several images obtained with HST/ACS. A closer inspection of those 
data reveals that it is actually a pair of stars, se parated by 0.1 5^^ Thes e stars were resolved 
in the F606W/F814W photometry performed by lAnderson et al.l (120081) . and placed i n thei r 
data base with F606W, F814W, V and / magnitudes. According to lAnderson et al.l (120081 ) . 
the brighter component of the pair (star ^^5273) has V = 17.170, while the fainter one 
(st ar #5274) reaches onl y V = 20.650. We note that the difference between V = 17.170 
of lAnderson et al.l (120081 ) and V555 = 17.37 of BPH originates from different photometric 
systems used by those authors. 

The HST archive contains four images taken with the F814W filter between 05:27 and 
08:04 UT on 2006-03-05 with Tg^p = 40 s (Program ID 10775, PI A. Sarajedini). From 
them, we extracted the photometry of both components of the pair and checked which one 
was variable. The star #5273 was constant to within observational errors (amounting to 
0.015 mag), while #5274 exhibited fluctuations with an amplitude of ~0.5 mag. The final 
proof that the fainter component of the pair is responsible for the variability detected in the 
ground-based data came from the analysis of a yet another set of 30 archival HST/ACS 
images, collected with the F435W filter on 2006-04-08 (Program ID 10615, PI S. Anderson; 
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all but three of these obtaine d with Texp = 340 s). We extracted the photometry with 
Dophot/Allstar (IStetson 1987), a nd tr ansformed it to the ACS- VEGA system, following 
the prescription of ISirianni et al.l ( 120051 ) . Again, the brighter star was remarkably constant 
(o"f435W = 0.013 mag), showing no periodicity with Pi or P2 (in fact, it showed no periodicity 
whatsoever). 



Based on V^- magnitudes given by lAnderson et al.l (120081 ) . we estimated the instrumental 



magnitude of the star #5274 in our ground-based reference image used for image subtraction. 
This allowed us to transform the ISIS light curve from differential counts to magnitudes 
(the curve was extracted at the position of star #5274, and the phase difference between 
ISIS template and ACS F606W image was accounted for). Next, we phased the V^-band 
magnitudes of the star #5274 using the ephemeris ([T]), averaged them using phase bins 
A(j) = 0.02, and interpolated them onto phases at which the F435W frames were taken. The 
interpolation enabled us to calculate the (V^ — F435W) color as a function of phase. We found 
that the maxima of the light curve were bluer by ~0.2 mag than the minima (see Fig. [1]). 
Such a behaviour is indicative of the irradiation effect, which produces a light curve with 
one maximum per orbital period, as is observed when Pi is used phasing. Alternatively, P2 
produces two maxima per period - a feature characteristic of the ellipsoidal effect (however, 
in this case the irradiation may also play a role). In Sect. [3] we argue that the actual orbital 
period Porb is most probably equal to Pi {^ 0.2628 d ~ 6.3 h). 



2.3. Chandra 

M4 was observed three times with the ACIS-S spectrometer on the board of Chandra 
X-ray Observatory on 2000 June 30 (Obs. ID 946), 2007 July 06 (Obs. ID 7447), and 
2007 September 18 (Obs. ID 7446) for 7.17, 12.63 and 13.31 h, respectively. During all 
observations the target was centered on the S3 chip, and the integration time was 3.24 sec 
per frame in the TIMED mode. All data were downloaded from the Chandra archive. The 
observations from 2000 June 30 were analyzed earlier by BPH, but we decided to repeat the 
analysis so that all data sets would be processed in exactly the same way. For the calibration 
of each observation the same version of CALDB 3.4.0 was employed. 

CXI was localized using the physical coordinates given by BPH (i.e. a = 16^23™34.12'^ 
and 6 = — 26°31'34.85"), and the light curves were exctracted using CIAO 4.3 software. All 
data sets were filtered for energies 0.5-7 keV, as suggested in CIAO ACIS analysis threads. 
Background flaring was accounted for with the help of CIAO tools dmextract and Ic-dean. 
Since the background was practically not flaring, good-time intervals did not differ much 
from total exposure times, amounting to 6.6, 12.06, and 12.96 h for observations 946, 7447, 
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Fig. 3. — X-ray light curves of CX 1. Left column: counts phased with the epehmeris 
([1]) without phase-folding, and binned in four intervals per period (filled triangles). Empty 
triangles: the same binning, but the bins are shifted by A0 = 0.125. Right column: counts 
phased with P2 used in eq. ([1]) instead of Pi, and binned in six intervals per period (filled 
triangles). Empty triangles: the same binning, but the bins are shifted by A0 = 0.083. No 
phase-folding was applied. Panels a, b and c show observations from 2000-06-30, 2007-07-06 
and 2007-09-18, respectively. = corresponds to the minimum of the optical light curve. 



and 7446, respectively. 

The source region was defined as a circle with a radius of 3" centered on CX 1, and 
the background was sampled in a surrounding annulus with inner and outer radius of 3"and 
10" (we checked that varying the outer radius within reasonable limits did not change the 
shape of the extracted light curves). Data sets 946, 7447, and 7446 yielded, 410±20, 610±24, 
and 494 ± 22 counts, respectively. The counts were binned in intervals of 0.25 and 0.167 
of the period for Pi and P2, respectively, and the resulting X-ray light curves are shown in 
Fig. O While rather noisy, they exhibit a modulation similar to that of the optical light 
curves, with a minimum near = 0. 

The same regions were used for the extraction of spectra with the specextract tool. To 
within standard errors, in each data set the same number of counts was detected as before. 
The spectra were binned, and each bin contained at least 15 counts. XSPEC v. 12.6.0 was 
used for fitting. We tried to fit three models: a power law, a black-body, and thermal 
bremsstrahlung, each convolved with a warm absorption model. For all three data sets the 
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Fig. 4. — X-ray spectrum of CX 1, derived from the combined exposures described in the 
text. The lower panel shows the residual to the power law fit. 



only reasonable fit was the power law. The black-body produced either rx^/dof > 2, or 
warm absorption dropping to zero, whereas thermal bremstralung required an unacceptably 
high temperature (199 keV). 

The power law fit to the 946 data set yielded an absorbing column density Nh = 
2.29 ±1.04 X 10^1 cm~2, and photon index T = 0.96 ±0.17 (xVdof = 25.5/22 = 1.16), which 
agrees well with the best-fit results of BPH, who obtained N^ = (1.95 — 2.86) x 10^^ cm~^ 
and r = 0.99 ±0.17 [x^/doi = 33.3/39 = 0.85). The absorbing column density derived from 
fitting also agrees with Nh = 2.18 x 10^^, converted from the visual ex tinction Ay = 1.2 2 



(Kalu zny et al. 2012, in preparation) accord ing to the p rescription of iPredehl fc Schmitt 
(119951 ). Assuming a distance to M4 of 2.2 kpc (JHarrislll996l ). we derived an unabsorbed X-ray 
luminosity i^o.5-2.5keV = 3.4 x 10^^ erg/s, and a hardness ratio -^0.5-1. 5kev/-^i.5-6.okeV = 0.21. 
Similar values were obtained for the remaining two observations. 

To improve the photon statistics, we combined all spectral data using the CIAO tool 
combine_spectra. The resulting cumulative spectrum of 1520 counts was binned in such a 
way that each bin contained at least 20 counts. The best-fitting model was again a power 
law convolved with warm absorption (Fig. H]), and its parameters did not differ much from 
those characterizing the individual spectra. We obtained Nh = 1-71 ± 0.4 x 10^^ cm~^, 
and r = 0.87 ±0.08 (xV^of = 67.84/61 = 1.11), with L0.5-2.5keV and L0.5-1.5kev/L1.5-6.0keV 
amounting to 2.66 x 10^^ erg/s and 0.19, respectively 

Our fit is significantly harder than that of BPH, resulting in a smaller hardness ratio. 
Their ratio is close to unity, but it was computed using correction factors averaged over 
three fits, whereas ours is derived from the best fit to the combined spectrum obtained from 



a much larger number of counts. 



Discussion and conclusions 



According to lZloczewskil (l2012l ). the star #5273 is a proper- motion member of M4. Upon 
averaging the subtracted ISIS-images on a seasonal basis, we found no evidence of bipolar 
residuals at the position of its weak companion. Such residuals are observe d for objects with 
notic eable proper motions with respect to the surrounding stellar field (JEyer fc Wozniak 
20011 ). Given the time base and the pixel scale of our ground-based images, we rule out 
a relative proper motion of the star #5274 in excess of 10 mas/y. Since on the proper- 
motion diagram m ost of the field stars are separated from the cluster members by 15 mas/y 
(jZloczewskil 120121 ) . this is consistent with the star #5274 being a member of M4. 



BPH tentatively classified CX 1 as a CV. However, between 1995 and 2008 average 
optical luminosity and amplitude of the star #5274 remained constant to better than 0.05 
mag. This would be very unusual for a CV whose substantial X-ray emission implies active 
accretion via LI overfiow. If CX 1 is indeed such a system, despite its stability, then its visible 
component should be a main-sequence star filling its Roche lobe, whose average density can 
be calculated from the formula 



p [g cm 
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107P-^ [hi 



(2) 



accurate to 3% for mass ratios 0.01 < g < 1.0 (lEggletonlll983l ). For Pi (^0.2628 d) and 
P2 (~0.5256 d) we get a p of 2.8 and 0.7 g cm~^, respectively. Since an 11 G yr old main- 
seque nce star with My = 7.15 and [Fe/H] = —0.7 should have p = 4.04 g cm~^ ( JDotter et al. 



20081 ). the orbital period of CX 1 cannot be equal to P2- -Pi might be considered marginally 
acceptable, were it not for the extraordinary stability of the star #5274. 

Another argument against CX 1 being a c ataclysmic variab l e follo ws from the ground- 



based UBV photometry of M4 published by iMochejska et al. 



(2002). Between the out- 



bursts, cataclysmic var iables typically have {U — V) ~ —1 (jWarnerlll995l ). In the data of 
Mocheiska etaP (j^OoJ the stars #5273 and #5274 are blended, and the {U - V) index of 
the unresolved blend amounts to 1.016, locating it right in the middle of the upper main 
sequence in the V/U — V diagram. With {U — V) = —1 (-0.6), the star #5274 would have 
an U of 19.55 (19.95) mag, and the color of the blend would be shifted S( U — V) = —0.2 8 



(-0.2 0). Such shifts would be easily seen on the V/{U — V) diagram of (IMochejska et al. 
2002h . 



We conclude that mean density, stability and color arguments taken together make the 
interpretation of CX 1 as a cataclysmic variable very unlikely. The only viable alternative is a 
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system whose primary is a neutron star instead of a white dwarf. If Porb = Pi then the optical 
variabihty of CX 1 must originate predominantly from irradiation (cf . Sect 12.21) . However, 
even in this case the star #5274 is likely to nearly fill its Roche lobe (first - it is heated, and 
therefore swollen; second - there must be a Roche-lobe overflow for the accretion to go on, and 
the X-rays to be produced). Alternatively, if Porb = -P2 then the ellipsoidal effect becomes 
the main cause of optical variations. In any case, the semi-detached configuration should 
be a reasonable approximation of the CX 1 system, and the estimates of the mean density 
of the star #5274 based on equation ([2]) remain valid. Having the same luminosity and 
temperature as similarly loc ated main-sequence stars of M4, the star #5274 must also have 
the same radius R ^O.GRq (JDotter et al.ll2008l ). For that size one gets a mass of ^OAMq 
and ~O.1M0, respectively for Pi and P2. Since the corresponding main-sequence mass is 
0.58 Mq, the Roche-lobe filling companion of CX 1 would have to be ablated either slightly 
(if Pnrb = -Pi) or very str ongly (if Porb = P2)- Simple experiments with the PHOEBE 
tool (jPrsa fc Zwitteiil2005l ) show that the V^-light curve for Pi can be easily recovered, e.g. 
in a system with neutron star of I.5M0, semimajor axis of 2.3P0, inclination of ~20° and 
illuminated-area temperature 20% higher than the rest of the photosphere (obviously, this 
does not prove that Pi is the actual orbital period, nor the semi-detached configuration is 
the only possible). 

For Porb = Pi there is only one X-ray minimum per period, which can be naturally 
explained as an eclipse, reminiscent of those observed in syst ems containin g millisecond 
pulsars J0024-7203W in 47 Tuc or J1740-5340A in NGC 6397 jHeinkel [201ll ) . Two X-ray 
mimima per period, appearing for Porb = P2 (see Fig. [3]), would be very hard to explain. 
We conclude that, while P2 cannot be ruled out, it seems much less likely than Pi. 

Thus, the available data favor the possibility that CX 1 is composed of a neutron 
star accompanied by a ^0.6Mq main-sequence star or a partially ablated, slightly oversized 
~O.4M0 star. If this interpretation is correct, then it must be a millisecond pulsar or a 
qLMXB. The X-ray luminosity and hardness ratio derived in Sect. 12.31 support the first 
possibility ; in fact, they l o cate CX 1 among the brightest millisecond pulsars shown in 
Fig. 6 of iBogdanov et al.l (120101 ) (the nondetection of the radio emission may be due to 
an unfavourable orientation). As for the radiation mechanism, the strongly nonthermal 
spectrum suggests that the dominant source of X-rays may be a relativistic shock, produced 



by th e pulsar wind interacting with the matter from the optical companion (IBogdanov et al. 



20051 ) . Future observations with XMM- Newton (more detailed X-ray light curves), HST 
(multicolor photometry) and VLT or Gemini South (IR spectroscopy with adaptive optics) 
should resolve the period ambiguity and more precisely establish the nature of this interesting 
object. 
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